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Abstract 



We derive a criterion for the breakdown of solutions to the Oldroyd- 
B model in M 3 in the limit of zero Reynolds number (creeping flow). 
If the initial stress field is in the Sobolev space iT m (R 3 ), m > 5/2, 
then either a unique solution exists within this space indefinitely, or, 
at the time where the solution breaks down, the time integral of the 
I/"") , L°°-norm of the stress tensor must diverge. This result is analogous to 

the celebrated Beale-Kato-Majda breakdown criterion for the inviscid 
Eluer equations of incompressible fluids. 

ON 
O 

o 

1 Introduction 

•T- 1 . 

X 

The Oldroyd-B model is a classical model for dilute solutions of polymers 
suspended in a viscous incompressible solvent [lj. Although it suffers, as a 
physical model derived from microscopic dynamics, from numerous shortcom- 
ings (e.g., polymers are allowed to stretch indefinitely), it is often considered 



■"Institute of Mathematics, The Hebrew University, Jerusalem 91904 Israel 
^Institute of Mathematics, The Hebrew University, Jerusalem 91904 Israel and CER- 
MICS, Ecole Nationale des Ponts et Chausses, 77455 Marne la Vallee, France 

■^Faculty of Mathematics and Computer Science, The Weizmann Institute of Science, 
Rehovot 76100 Israel and Department of Mathematics and Department of Mechanical and 
Aerospace Engineering, University of California, Irvine CA 92697-3875 



1 



BEALE-KATO-MAJDA CRITERION FOR AN OLDROYD-B FLUID 2 



as a prototypical model for viscoelastic fluids, and has therefore been the 
focus of both analytical and numerical studies. 

At present, there is no global- in-time existence theory for the Oldroyd-B 
model. The notable difference between the Oldroyd-B model and its Newto- 
nian counterpart, the incompressible Navier-Stokes equations, is that in the 
viscoelastic case, global-in-time existence has not even been established in 
two dimensions nor in the creeping flow regime, i.e., when the momentum 
equation is a Stokes system. The reason for this difference can be understood 
by observing structural similarities between the inertia-less Oldroyd-B model 
and the Euler equations (in three dimensions), or the 2D quasi-geostrophic 
flow equations (in two dimensions) [2]. 

Since the early 1970s, numerical simulations of the Oldroyd-B model (as 
well as other viscoelastic models) have been infested by stability and accuracy 
problems that arise at frustratingly low values of the elasticity parameter 
(the Weissenberg number) [21 H]. While some of these difficulties have been 
elucidated [5], it is to a large extent still a mystery why computations break 
down in the low-Reynolds-high- Weissenberg regime. As is often the case in 
such situations, numerical data are by themselves not sufficient to explain 
the reasons for this breakdown. In the absence of a well-posedness theory, 
it is not even clear in what spaces solutions have to be sought. Thus, the 
development of such a theory is of major importance for both theoretical and 
practical purposes. 

This situation is analogous to that of incompressible Newtonian fluids 
at high Reynolds number, where global-in-time existence has not yet been 
established. For a Newtonian fluid (in three dimensions), however, there is 
a prominent observation due to Beale, Kato and Majda (BKM) [6], which 
states a necessary and sufficient condition for the breakup of solutions at 
finite time. Specifically, the Euler equations of incompressible inviscid fluids 
in vorticity formulation take the form 

^ + (« • V)w = (Vu) V (1.1) 

with initial condition cj(-, 0) = ujq. Here oj = V x u is the vorticity and u is 
the velocity field. The BKM theorem states that if ujq belongs to the Sobolev 
space H m {R 3 ), m > §, then either there exists a solution u>(-, t) e H m (M 3 ) for 
all times, or, if T* is the maximal time of existence of a solution in if m (R 3 ), 
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then 



and in particular, 



lim / \\u>(-, s) || r,°° ds — oo, 

t/T* 



lim sup ||u>(-, £)||l°° = oo. 

t/T* 



That is, the breakup of solutions in any Sobolev norm necessitates the di- 
vergence of the L°°-norm of the vorticity. The practical implication of this 
theorem is that breakdown cannot be attributed, say, to the failure of some 
high derivative. The blowup of the vorticity itself, in the supremum norm, 
is the signature of any finite-time breakdown. For another criterion of sin- 
gularity formation see [TJ; for up-to-date surveys see [8j[9]. 

The goal of this paper is to establish a similar result for the three- 
dimensional Oldroyd-B model, in the zero-Reynolds number regime. In 
this regime, a closed equation can be written for the polymeric stress field 
cr = <t(x, t); this equation is similar to the vorticity equation (11. II) . We start 
by establishing the local-in-time existence of solutions in any Sobolev space 
if m (R 3 ), m > 5/2. Following then the approach of BKM, we prove that if 
the initial stress is in H m (W. 3 ), then either a solution exists for all time, or, 
if T* is the maximum existence time, then 



lim / ||<t(-, s)||l«> ds = oo. 
t /" T * Jo 

This result is independent of the Weissenberg number, and in fact, holds even 
if one sets the Weissenberg number to be infinite. From a theoretical point 
of view, this breakdown condition implies that global-in-time well-posedness 
hinges on an a-priori bound for the supremum norm of the stress. 

Recent work along these lines comprises a BKM-type analysis by Chemin 
and Masmoudi [10J. The notable difference between their analysis and the 
present work is that they treat the Oldroyd-B model including inertia, how- 
ever, their analysis is restricted to two-dimensional flows. Lin et al. [TTj 
analyze the inertial Oldroyd-B model without the relaxation term (infinite 
Weissenberg number) and establish global-in-time existence for small ini- 
tial data. Finally, Constantin [12] studies a class of kinetic models in the 
form of a Stokes equation coupled to a nonlinear Fokker-Planck equation, for 
which he proves global-in-time existence. The extension to a two-dimensional 
Navier-Stokes system coupled to a nonlinear Fokker-Planck equation appears 
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in Constantin et al. [13] and Constantin and Masmoudi [14]. In both cases 
the analysis benefits from the fact that the polymeric stress remains bounded 
by construction, which as our analysis shows, is the key to well-posedness. 



2 The model 

The Oldroyd-B model describes a fluid in which polymer molecules are sus- 
pended in a viscous incompressible solvent. The equations of motion in the 
creeping flow regime are 

= — Vp + u s Au + div cr 

diva = (2 1) 

^ + {u ■ V)a - {Vufa - cr{Vu) = -\a + ^ [Vu + (Vuf] , 

ut A A 

where u is the velocity of the fluid, p is the hydrostatic pressure, cr is the 
extra-stress tensor due to the polymer molecules, v s is the solvent viscosity, 
v p is the polymeric viscosity, and A is the elastic relaxation time. The velocity 
gradient is defined with components (Vu)y = duj/dxi. 

The first two equations in (12 .ip are a Stokes system for an incompressible 
fluid, whereas the third equation is the Maxwell constitutive equation for the 
extra-stress pQ. The flow is assumed to take place in the unbounded three- 
dimensional space M 3 . Initial data need only be prescribed for the stress, 
tr(x,0) = (T (x). 

The system (12 .ip can be turned into a closed equation for cr, by solving 
the Stokes system, and expressing the flow field in terms of the stress field. 
Specifically, the induced velocity field is given by 

U A X ) = I M u ] {y)dk<?ki{x - y) dy, 



where 



M$\y) = -^ + yjVl 



\y\ w 



is the Stokes kernel (see Galdi [15] . pp 189-195). Using integration by parts, 
it can be rewritten as 

Uj (x) = J- / MjH(y)a k i(x - y) dy, (2.2) 

OTlfg J R 3 
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where 

M jkl(y) = --^r + ^r- 

Here and below we adopt the Einstein summation convention, whereby re- 
peated indexes imply a summation unless otherwise specified. In tensor no- 
tation we write (12.2ft as 

u{x) = J- / M«(y) : <r(a: - y) dy, (2.3) 

where for 2-tensors a, b, the : product is defined by a : b = tr (a T b). We 
then rewrite the constitutive equation as a closed evolution equation for the 
stress field, 

^ = F(<t), <t(-,0) = <t , (2.4) 

where 

F{tr) = -(u ■ V)<r + (Vu) T er + <r(Vu) - t ct + "V [ Vw + ( Vn )T » ( 2 - 5 ) 

A A 

and u is given by ( 12.31) . Equation ( 12.41) is viewed as an evolution equation or 
an ordinary differential equation (ODE) in an infinite-dimensional function 
space. We observe that the solution cr is a symmetric tensor whenever <x is 
symmetric. 

For later use, we derive the linear relation between the stress field cr and 
the velocity gradient Vu, obtained by differentiating (12. 3p and integrating 
by parts. This yields a singular integral (the integrand is a homogeneous 
function of degree —3 that averages to zero on the unit sphere), from which 
one has to extract the singular part, 

d iUj {x) = -— (cxijix) - -8 ij tTa(x))+-— (P. V. ) / M[f kl {y)a kl {x-y) dy, 
bv s \ 3 J mu s J R3 

where 

,,(2)/ \ SijSki ViVjhi + 2y j y l 6 ki + kjVkVi . l^yiyjykyi 

\y\ \y\ \y\ 

and (P. V. ) stands for the principal value of a singular integral. In tensor 
notation, 

V«(x) = --L(o-(aj)-|trff(x))+-^-(P.V.) / M®(y) : cr(x-y) dy. 

bu s \ 3 J 8irv s J R 3 

(2.6) 
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3 Local-in-time existence 



In this section we prove the local existence and uniqueness of solutions to 
the Oldroyd-B model (12.11) . The proof is based on energy methods, and 
closely follows the existence proof for the Euler and Navier-Stokes equations 
in Majda and Bertozzi [IB] . Differences between the two cases are highlighted 
along the treatment. 

We denote by H m (R 3 ) the Sobolev spaces of scalar, vector and tensor 
fields in R 3 . The corresponding norms are denoted by || • || m and are defined 
by 




where here / denotes either a scalar, a vector or a tensor field and a = 
(ai,«2,ct3) is a multi-index of derivatives; the norm || • || is the L 2 (IR 3 )- 
norm. For example, if / = / is a tensor field with components then 



\D"f\\l 



0H 



dx" 1 dx% 2 dx 



a 3 A? 1 



X 



dx. 



We denote by || • and || • H^oo the L 9 (IR 3 ) and L°°(R 3 ) norms, respectively. 
Throughout this section and the next one, we use the symbols C, K, to denote 
either positive constants, or, depending on the context, bounded functions 
of their arguments. 

Our local-in-time existence theorem is: 



Theorem 3.1 Let cr G H' m (R 3 ) for m > 5/2. Then there exists a time 
T > depending on ||cr || m only, so that the Oldroyd-B system (12.11) . or 
equivalently the Hilbert space-valued ODE ( 12.41) has a solution <r in the class 

ct e C([0,T];iJ m (M 3 ))nC 1 ([0,T];i/ m - 1 (M 3 )). (3.1) 

Observe that system (12 .ip is time- reversible, hence Theorem 13. II is also valid 
backward in time. 

Since the proof is long and technical, we describe here its outline, and 
prove each step in a separate subsection. Standard definitions and inequali- 
ties are grouped in Appendix [XI 
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Subsection 13.11 We start by constructing smooth approximations to a. 
We consider a mollified version of (12.41) . 

^£ = F £ (o- e ), rr e (.,0) = <r„, (3.2) 
where e > is the mollification parameter, 

F £ ((T £ ) = -J £ [U £ ■ V(J £ <T £ )] + (VU £ ) T (T £ + (T £ V(U £ ) 
~ \*e + ^ [VU £ + ( Vu E f] , 

u £ is given by the integral (12.31) with a replaced by cr £ , and the mollification 
operator J £ is defined by (1A.1I) . We then prove 

Proposition 3.2 Let <x e if m (M 3 ) for m > 3/2. Then there exists a time 
T £ > depending on ||<To|| m only, so that (13. 2p /ias a unique solution 

<T e eC\[0,T e y,H m (R 3 )). 

The main reason for introducing the mollified equation is to enable us to 
use classical theory of evolution equations (ODEs) in Banach spaces to prove 
short time existence of unique solutions. Most importantly, the solutions of 
the mollified equation are regular enough, which is enabling us to use classical 
tools for deriving estimates without a need for additional justification. 

Finally, we also observe that by the uniqueness of the solution of (13. 2)) , 
the tensor <t £ is symmetric whenever cr is symmetric. We will be using this 
fact later in our estimates. 

Subsection 13.21 Using the continuation theorem for autonomous ODEs, 
and an a-priori estimate of the form 

sup ||cr £ || m < if(||er || m ,T), (3.3) 

0<t<T 

where T = T(||cr || m ) is independent of e, we show that the family of mollified 
solutions <j £ can be continued uniformly up to a common time T. Here we 
need a slightly higher degree of regularity: 

Proposition 3.3 Let cr G H m (M. 3 ) for m > 5/2. Then there exists a time 
T = T(||<7o|| m ) > independent of e, such that the mollified equation (13.21) 
has a unique solution cr £ G C 1 ([0, T], H m (R 3 )), satisfying the uniform bound 

(ED. 
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Subsection 13.31 We show that the family of mollified solutions, cr £ , forms a 
Cauchy sequence in C([0, T], L 2 (R 3 )), hence strongly converges to a function 
TeC([0,T],L 2 (I 3 )). 

Proposition 3.4 Let ctq G H m {R 3 ) for m > 5/2. Then the family of molli- 
fied solutions cr £ G C 1 ([0, T], H m (M 3 )) forms, as e — > 0, a Cauchy sequence 
in C([0, T], L 2 (M 3 )), /ience converges to a function, which we denote by cr. 
Moreover, for every < t < T , we have cr(-,t) G H m (R 3 ), and cr satisfies 
the same bound 

sup ||cr|| m < if(||<x || m ,T), 

0<t<T 

as the family of mollified solutions. 



Subsection 13.41 Using the technique of interpolation we show that a e 
strongly converges to cr in all intermediate norms C([0, T], H m (M 3 )), < 
m' < m. We then proceed to prove continuity of the limit in the highest 
norm, cr G C([0, T], H m (R 3 )). 

Proposition 3.5 Let tr G H m (R 3 ) for m > 5/2. Then the C([0, T], H m ' (R 3 )) 
limit cr of cr £ , for every m! G (0, m) is continuous in cr G C([0, T], H m (M 3 )). 



Subsection 13.51 We finally show that cr, the limit of cr E , is a solution of 
(12.41) in the space 

cr G C([0,T},H m (R 3 ))nC 1 ([0,T},H m - 1 (R 3 )). 



3.1 Local-in-time existence of mollified solutions 

In this subsection we prove Proposition 13.21 We approximate the Oldroyd- 
B system (12. ip . or equivalently the Hilbert space- valued ODE (evolution 
equation) (12. 4p by a mollified equation for a smooth approximation cr e of cr, 



dt 
where 



<la - F £ (cr £ ), cr £ (;0)=cr , (3.4) 



F £ {cr £ ) = -J £ {u £ ■ V(J £ cr £ )] + (Vu £ fcr £ + cr £ (Vu £ ) 

1 v v r _ , T1 (3-5) 

- jcr £ + [Vu £ + (Vu £ ) T ] , 
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u £ is given by 112. 3ft with cr replaced by cr e , and the mollification operator 
J £ is defined by 1 1A.1I) in the appendix. Comparing with 112.51) . we note 
that mollification is only used in the advection term. As will be shown, 
the gradient of u e has the same degree of regularity, with respect to the 
if m (R 3 )-norms, as cr e , hence no additional mollification is needed. 

To prove that 113 ,4p has a local-in-time solution we use Picard's theorem 
over Banach spaces. Specifically, we work within the Banach space H m (M. 3 ) 
with m > 3/2. Picard's theorem for functional evolution differential equa- 
tions states that if there exists an open subset O C H m (M. 3 ) such that 

1. F £ : O -> H m (R 3 ). 

2. F £ is locally Lipschitz continuous, i.e., for any cr e O there exists an 
open neighborhood of cr, U C O, and a constant L > such that for 
every T\, t 2 G Z7 

||F e (ti) - F e (r 2 )|| m < L||ri - r 2 || m 

then there exists, for every cr 6 O, a time T e > and a unique solution 
cr e G C 1 ([0,T £ );O) of (E3J). 

Two properties that are being used extensively throughout this section 
are: (i) the Calderon-Zygmund (CZ) inequality flA.lip (see Appendix), from 
which it follows that cr £ G H m (R 3 ) implies that the velocity gradient Vu E , 
and a-fortiori the velocity u £ itself, are in H m (M. 3 ) as well, 

||«e||m < — \\Ce\\m , ||Vt* e || m < — ||cr £ || m . (3.6) 

(ii) For m > 3/2, H m (M. 3 ) forms a Banach algebra, i.e., 

\\fg\\m < c||/||m||fl , ||m- 

Combining these two properties with the smoothing properties ( 1A.3}) — 
( 1A.4P of J £ (see Appendix), it follows at once that F £ is a mapping if m (M 3 ) — > 
H m (R 3 ). We set 

O = {* E H m (R 3 ) : H| m <r}, 

where r is sufficiently large such that cr G O. It remains to show that there 
exists a positive constant L = L(r), such that F £ is Lipschitz continuous in 
O with constant L. 
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To avoid lengthy expressions, we split F £ into a sum of four terms, 



F 1 + F 2 + F 3 + F 4 , 



where 
Fi(tr e 
F 3 {(T £ 



J £ [u £ ■ V(J £ a £ )] 



X 



F 2 ((T £ ) = (VU £ ) T (T £ + (T £ (VU £ 

FMs)= ^[Vu £ +(Vu £ f] 



and show that each of these four terms is Lipschitz continuous. That is, let 
Ti,T2 £ O and let u x , u 2 be their corresponding velocity fields satisfying the 
Stokes system ( 12. 31) ; we show that each of the Fj verifies a bound of the type 



|F 7 -(ri)-F 7 -(r 2 )|L<L(£, llnlL 



T 2 



Tl - T 2 



where L is a monotonic function of its last two arguments, and hence it is 
bounded by L(e,r,r). 

The Lipschitz continuity of the linear function F% is trivial. The Lipschitz 
continuity of F± follows from the CZ inequality (lA.llj) . which implies 



\Vu 2 - V«i| 



c „ 

<- r 2 



Til 



For the advection term F\ we have 

\F x {t 2 ) - F x {r x )\\ m = \\J £ [u 2 ■ V(J £ r 2 )] - J e [u x ■ V{J e r x )\\\ m 

< C \\u 2 ■ V(J £ r 2 ) - u x ■ V(J £ n)|| m 

< C \\u 2 ■ V(J e (r 2 - n))|| m + C \\(u 2 - u x ) ■ V{J £ r x ) 



where in the passage from the first to the second line we used (IA.3j) with 
k — 0, and in the passage from the second to the third line we added and 
subtracted equal terms and used the triangle inequality. Using the Sobolev 
calculus inequality (IA.6I) . 



\F x (r 2 )-F x (r x )\\ m <C 



w 2|| L oo ||V(J £ (r 2 - T X ))\\ m 

+ ||V(J £ (T 2 - T x ))\\ LOC \\u 2 \\ m + \\U 2 - Ui|| £o 
+ ||V(J e Ti)|| ioo \\U 2 - Will, 



iVfXr 



ei 1 m 
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Using then property (1A.3I) of the mollification operator and the Sobolev em- 
bedding (1A.9jl . we obtain: 

C 

||F 1 (r 2 ) - -Fi(ri)|| m < — (||it 2 || m ||r 2 - Ti|| m + ||ri|| m ||w 2 - u x \\ m ) . 
Finally, using CZ inequality 

||F!(t 2 ) - Fi(Ti)|| m < C(e, ||Ti|| m , ||T 2 || m )||T 2 - Ti|| m . 
Remains the deformation term F 2 . Using the Banach algebra property flA.5j) . 

||F 2 (r 2 ) - F 2 (n)|| m = ||(Vtii) r T 1 + Ti(V«i) - {Vu 2 ) T r 2 - r 2 (Vw 2 )|| m 

< 2||n(Vwi) -T 2 (Vu 2 )\\ m 

< 2||Vui|| m ||r 2 - ri|| m + 2||r 2 || m ||Vu 2 - Vni|| m . 

One more application of the CZ inequality gives, 

||F 2 (t 2 ) - F 2 (r 1 )\\ m < C(||ri|| m + ||r 2 || m ) ||r 2 - Ti|| m . 

This shows that F £ is locally Lipschitz in O, hence local existence for the 
mollified equation (13.41) . 

3.2 Energy estimates for the mollified solutions 

In this subsection we prove Proposition 13.31 i.e., that the mollified solutions 
can be continued uniformly up to a time T > that does not depend on 
e. To do so, we first obtain an a-priori estimate, whereby if the solution a £ 
exists up to time T, then 

sup \\cr £ (-, t)\\ m < C(T, \\cr \\ m ) = K, 

0<t<T 

where all the <t £ have the same initial condition. The existence of the solution 
up to that time follows from the continuation theorem for autonomous ODEs 
in Banach spaces. Taking for domain 

0= {(re H m (M. 3 ) : ||cr|| m < K) , 

the solution either exists up to time T, or leaves the set O before that time, 
which we would have ruled out by the above estimate. 
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To obtain an a-priori uniform bound on ||<r £ (-, t)\\ m , we use an energy es- 
timate. Starting from the mollified equation ( 13.41) . we take its a-th derivative 
(\a\ < m), and then an inner product with D a cr £ . This yields an "energy" 
equation, 

~±\\D<*<T E g + j\\D a cr £ \\l = - (D a er £ , D a J £ [u £ ■ V(J £ a £ )}) 

+ 2 (D a er £ , D a [a £ (Vu £ )]) + ^ (D a er £ , D a (Vu £ )) ^ 



h + h 



'3, 



where we used the symmetry of a £ in the last two terms. 

We observe that cr £ is as smooth as the initial data cr 6 if m (R 3 ), and 
J £ cr £ G C°°(IR 3 ), therefore the above estimates should be interpreted in the 
strong sense. In particular, the time derivative is classical and we do not 
require additional justification. 

Since we need a bound that does not depend on e, we cannot use the 
smoothing properties of J £ . On the other hand, we are not concerned by 
finite-time blow up as long as the time horizon is independent of e. 

We start with the advection term. Using the fact that J £ commutes with 
weak derivatives and is symmetric, we write 

-h = (D a (J £ a £ ), D a [u £ ■ V(J £ a £ )}) . 

We then add and subtract from the second argument of the inner product 

u £ -D a V(J e <r E ), 

whose inner product with D a (J £ cr £ ) vanishes due to the incompressibility of 
the flow, i.e., due to u £ being divergence-free. Using the Cauchy-Schwarz 
inequality and (1A.3j) for k = we obtain 



h < C||D £ V e || ||D a K • V(J £ <r £ )] - u £ ■ D a V(J £ (T £ )\\ . 

We then invoke the Sobolev calculus inequality (IA.7I) to get 

h < C||£> Q <r £ || [||Vu e ||£~||V(J £ <r £ )|| m _i + K|| m || V(J £ <r £ )||^] . 

Using repeatedly the property flA.3j) of J £ , the Sobolev embedding theorem 
f)A.9j) . and the CZ inequality, we get the bound 
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Note that in order to bound ||V<T e ||^oo by ||cr £ || m we need m > 5/2. 

We turn to I2, where we use the Cauchy-Schwarz inequality, the Banach 
algebra property of H k (R 3 ) for k > 3/2 and the CZ inequality, 



I 2 = 2 (D a tr e ,D a [tr e (Vu E )]) < 2||<r E || m ||<r 6 (Vti. 



e J \\m 



„ ,.9., Cn ,.<t (3-9) 

<CKOV« e |U<-||<r e ||^. 

Remains J 3 , which we estimate using the Cauchy-Schwarz inequality followed 
by the CZ inequality, 

h = ^ (D a cr £l D a Vu £ ) < p^Holl^V^llo < C^\\*sf m . 

(3.10) 

Combining the three estimates (13.81) . (I3.9P and (13.101) we obtain an energy 
inequality 



d ., 1 Up C2 ,|2 

I Fe m + T Ff m — °e m H CT e m 

at A Az/o z/o 



from which we conclude the existence of a time T = T(||cr || m ), independent 
of e, for which all the er £ exist and have a common bound in H m (M. 3 ), 

sup ||<T e || m < rr ^ i ° T " n = K(||o- || m ,T), (3.11) 



o<t<r c 3 + c 4 (l - e csl ) cr 



where c 3 = A 1 {c 1 v p /u 3 - 1), c 4 = c 2 /z/ s and T < j^jlog(l + |c 3 |/c 4 ||cr || m ). 
This concludes the proof of Proposition 13.31 

Remark: Substituting back the dimensional parameters, our expression 
for the uniform existence time is 

A u s \ciu p /u s - 1| 



T< i — / — n lo § 1+ mi 

\civ p /v s -l\ \ c 2 X\\cr 

First, it follows that this time remains finite in the A — > 00 limit (i.e., infinite 
Weissenberg number), in which case T < z^/c2||cr || m . Second, this time is 
unbounded (hence, global-in-time existence follows) if v s jv v > c 2 and the 
initial data are sufficiently small, namely ||cr || m < |c 3 |/c 4 . 
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3.3 Convergence of <r £ in C([0, T], L 2 (R 3 )) 



We proceed to prove Proposition I3.4[ whereby the sequence cr £ forms, as 
e — > 0, a Cauchy sequence in the space C([0, T]; L 2 (R 3 )). Here, T is the 
uniform existence time established in the previous subsection. Specifically, 
we show that for cr £ ,cr £ t £ if m (IR 3 ) solutions of (13.41) with the same initial 
condition cr , the following holds: 

sup ||<r E - <T e /|| < C(\\er \\ m ,T)max(e,e'). 

0<t<T 

Hence follows the existence of cr £ C([0, T], L 2 (M 3 )), such that 
lim sup \\(T £ (-,t) - tr(;t)\\ = 0. 



0<t<T 



As in the previous subsection, we start with an energy equation, this time 
for the difference er £ — (T e r. 



1 d II 112 , 1 || in 

2^K-^|lo + x K-^llo 

- (J £ [u £ ■ V(J £ cr £ )] - J £ ,[u £ , ■ V(J E /0V)], cr £ - cr £ , 



(3.12) 

+ 2 (cr e Vit e - cr e / Vu E >, a £ - oy) 

+ (Vu £ - Vw £ /, <r £ - oy) = Ii + h + J 3 - 

^2 and I3 are easily estimated by the same manipulations as in the pre- 
vious subsection, 

h < 2||cr e Vw e - t7yVw £ /|| Q ||<x e - oy|| 

= 2||(<x £ - oy)Vu £ - a £ ,(Vu £ , - Vu £ ) || ||cr £ - oy || 

< 2 [\\cr £ - 0v||o||Vu e ||i,« + ||cr £ /|| L oc||Vn £ / - Vu e || ] ||cr £ - oy|| 

< C(||Vtt e ||ioo + ||o- £ /||ioo) \\cr £ - cr £ ,|| 2 

< C(||c"e||m + ||c£'||m)||cr£ — OV||o> 

(3.13) 

and similarly, 

h < ^||Vw £ - Vw £ /|| ||cr £ - oy|| < C\\a £ - ovU 2 . (3.14) 
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Remains the advection term I\, which we first split as follows, 

-h = ((J £ - J £ ')[u £ ■ V(J £ cr £ )),cr £ - a £ >) + {J £ >[{u £ - u £ >) ■ V(J £ <r £ )], a £ - c 
+ (J £ >[u e , ■ V(J £ - J £ >)a £ ),a £ - oy) + (J £ ,[u £ ■ VJ e /(cr e - oy)],<x e - <x £ , 

The last term vanishes because u £ is divergence-free. For the first three terms 
we use the Cauchy-Schwarz inequality, obtaining thus 

71 < \\{Je - Je>)[u £ ■ V(J £ (T £ )}\\o + || J £ >[(u £ - U £ >) ■ VJ £ a £ }\\ 

We - O"e'||0 

+ \\J £l [u £ , ■ V(J £ - J e 0<r e ]||o = Ax + A 2 + A 3 . 

(3.15) 

By (TOl) the outer J £ , can be replaced in A2, A3 by a constant prefactor. 
^2 is easily estimated by 

A 2 < C\\u £ - ii e /|| ||V(J E cr £ )|| L oo < C\\er £ - CT £ /|| ||<r e || m , (3.16) 

where we have used the CZ inequality and the Sobolev embedding theorem, 
with m > 5/2 (see Appendix). 

Note that the two remaining terms have the factor (J £ — J e /), which is 
"small" in the following sense. By flA.2h follows that 

\\(Je ~ Je>)f\\ < \\(Je ~ /)/||o + We' ~ I)f\\ < C||/||l max( £)£ '). 

Thus, Ax can be estimated by 

Ax < C\\u £ ■ V(J £ <r £ )||i max(£, e') < C||<t £ ||^ max(e, e'), (3.17) 

where the last inequality follows from the very rough estimate of the H 1 - 
norm by the H m ~ l norm, and the Banach algebra property of H k (M. 3 ) for 
k > 3/2. 

y4 3 verifies an estimate similar to Ax- Gathering the expressions for Ax, 
A 2 , A 3 , J 2 , I3, and substituting into the energy equation (13.121) . 

Tj-^ll^-o'e'llo < C , ||cr e ||^max(e,eOlk £ -CT £ nio+C'(||cr e ,|| m +||cr e || m )||cr £ -cr e ,|| 

(3.18) 

We now use the uniform bound A3. 11 j) to obtain 

^||cr £ - ov||o < C(||cr || m ) [max(e,e') + ||<r e - cr e /|| ] , 
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which upon integrating yields 

sup ||<r e - ov ||o < e c{ll(Tollm)T max{e,e'), (3.19) 

0<t<T 

and we used here the fact that er e and ov satisfy the same initial conditions. 
Therefore, cr e is a Cauchy sequence in the Banach space C([0, T}; L 2 (R 3 )) 
and hence it has a limit cr G C([0, T]; L 2 (IR 3 )). In particular, (13.191) implies 
that for the limit er we have 

sup ||0 e - 0|| o < e c(ll<TolU)T £. 

0<t<T 

The uniform boundedness (13.111) of the cr e implies by the Banach-Alaoglu 
theorem that for every t <T the sequence cr £ (-,t) has a subsequence that 
converges weakly in H m (M. 3 ). This limit must however coincide with the L 2 
limit, er(-,t). Moreover, the Banach-Alaoglu theorem also implies that 

\\cr(-,t)\\ m < Hminf ||0 £ (-,t)|| m , 

£—►0 

thus it follows from (13.111) that for every t G [0, T], 

sup \\er{-,t)\\ m <K{\\cr \\ m ,T). (3.20) 

0<t<T 

Note, however that we do not yet know that er is a continuous function from 
[0, T] into if m (R 3 ). It is the task of the next subsection to show that cr is in 
C([0,T},H m (R 3 )). 



3.4 Continuity in H m (R 3 ) 

In this section we prove Proposition 13. 51 whereby a G C([0, T], H m (M. 3 )). We 
do it in several steps. 

We start by showing that the mollified solutions u e converge to cr in all 
"intermediate" norms, C([0, T], H m '(R 3 )), for all < m' < m. For that we 
invoke the following interpolation lemma in Sobolev space: 

ll-T-ll <r r ll-i-ll 1 ~ m '/ m ll->-ll m 7 m 

||T||m' — ' || II ' II m ) 

valid for all r G H m (M. 3 ) and < m' < m. Substituting r = cr £ — cr, 

II /r /-rll <- n Wrr _ 1 1 1 / m 1 1 ^ ||m'/m 

H^s — ^Wm' S O m ||0 e — 0|| o ||CT e — CT|| m ' . 
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Using the uniform boundedness (13.111) of cr £ and (13.201) . 

K-<r||™' /m < (KIL + \\<T\\m) m ' /m < [2if(||cro|| rai r)]™v™ 

thus we obtain 

sup K(-,t)-<r(-,t)|U <{2K{\\<T Q \\ m) T)] m 'l m sup \\cr £ {;t)-cr{;t)\\l~ m ' h 

0<t<T 0<t<T 

i.e., uniform convergence er e —>■ cr in all intermediate norms. 

To show that cr is time- continuous in the highest norm, we first show that 
cr is time-continuous in the weak topology of H m (M. 3 ). That is, we claim that 
for every G H- m (R 3 ), 

g{t) = (0,<x(-,t)> 

is continuous in time, where (•, •) is the dual pairing between H m (M. 3 ) and 
H~ m (R 3 ). Since cr is time-continuous in the strong topologies of all the in- 
termediate norms, it follows that (<p, cr(-,t)) is time- continuous for all <fi G 
H- m '(R 3 ), but since the latter is dense in H m (M 3 ) and cr satisfies the uni- 
form bound (13.201) in H m , the continuity of g(t) follows. 

As is well-known continuity in the weak topology of a Hilbert space sup- 
plemented by the continuity of the norm yields continuity in the strong topol- 
ogy. Thus, it remains to show that ||cr(-,t)|| m is time-continuous. 

We start by showing continuity at the initial time t = 0. For ft > we 
have 

\\tr(-, h) - <7o||m = \\<r{; h)\\ 2 m - \\cr \\ 2 m - 2 (cr(-, h) - cr , cr ) m 

where (•, -) m denotes the inner-product in H m (M 3 ). As h — ► + , the last term 
vanishes by the time-continuity of cr in the weak topology in H m (R 3 ). This 
yields, 

\\cr \\ m < limmf \\cr(-, h)\\ m . (3.21) 

To obtain the reverse inequality we observe that by following similar steps 
as in the above proof on can obtain the following modification of (13.201) . 

sup ||<r(-,t)|| m < K(\\cr \\ m ,T), 

0<t<T 

for all r G (0, T], where K is given by (13.201) . Therefore, by taking r = h 
and assuming h is small enough, the above inequality and (13.201) yields, 

||o"(-> h)\\ m < ||o"o|| OT + Ch \\cr \\ m (1 + ||cr || m ) . 
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Since ||er(-, h)\\ m is bounded we may let h — ► + , obtaining 

limsup || cr(- , h)\\ m < ||cr || m , 

h\0 

which together with (I3.2ip implies right-continuity at t = 0; left-continuity at 
t = follows from the fact that the Oldroyd-B system can be time reversed 
(unlike parabolic equations such as the viscous Navier-Stokes equations). 

It remains to show that a is time continuous at any arbitrary time s G 
[0, T}. We use the fact that cr(-,s) G H m (R 3 ) to construct a new set of 
mollified solutions with initial data at t = s, <x E (-, s) = <x(-, s). By the same 
line of reasoning as before, these solutions converge in C([0, T], L 2 (M 3 )) to a 
solution cr which belongs, at all times, to H m (R 3 ) and is continuous at the 
initial time t = s. Let cr £ (-,t) be as before. Now, one can follow the same 
steps as in Section 13.31 to show that 

2^\\ a s - &e\\l < C (\\a £ \\ m + \\a £ \\ m )\\er £ \\ m £\\a £ - d- £ \\ 

+C(\\<7e\\m + ||cr e || m )||cr e — o- e ||o 
(cf. (I3.18P ). Similar steps lead that around t = s 

W&eWm < K(\\cr £ (-,s)\\ m ). 

Therefore, we have 

^||cr £ - cr £ || < C(||o- e (-,s)|| m ) (e+ \\tr e - & e \\ ) . 
Integrating, we obtain 

||<r e (., t)-a £ (; t) ||o < e C(ll^(-.-)IU)(*-) H^., S )-& £ (; S) || 0+e C(H^(-.«)IU)(*-^) £ . 

Now, as we let e — > 0, we know that er £ (-,s) — > cr(-,s), and we know that 
cr e (-, s) = cr(-, s), thus we conclude 

||<T(.,*)-cr(.,*)|| o = 0. 

Since <r is continuous at t = s so is cr. 
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3.5 <x is a solution of (12.41) 

Having shown that cr G C([0, T], H m (M. 3 )), it remains to show that cr is 
indeed a solution of (12.41) . and that 

cr E C([0,T],H m (R 3 ))nC 1 ([0,T},H m - 1 (W i )). 

To show that, we refer once more to the mollified solutions, whose evolu- 
tion satifies the integral equation, 



t) 



r + f F £ (cr E (-,s))ds. 
Jo 



We now exploit the convergence of cr £ to cr in all the intermediate norms. 
Specifically, we set 5/2 < m! < m, and claim that 

<r e -<r mC([0,T},H m '(R 3 )) 

implies that 

J £ [u £ ■ V(J £ cr £ )] -> u ■ Vcr 
Vm £ -> Vtt 

in C([0,T],fl' m '- 1 (R 3 )). The last two identities follow from the CZ inequal- 
ity (lA.lip and the Banach algebra property of H k (R 3 ) for k > 3/2. The 
convergence of the advection term follows from the same considerations, up 
to the loss of one order of regularity due to the gradient of <r £ . Thus, 



cr(-,t) =cr + [ F(tr(;s))ds, 

Jo 



which proves that cr £ C 1 ([0,T],if m 1 (M 3 )), and satisfies the differential 
equation (|2.4p . 

Remark: Based on the previous remark, this also proves global-in-time 
existence for small initial data. 



4 A Beale-Kato-Majda breakdown condition 

Having proved the local-in-time existence of solutions to the Oldroyd-B equa- 
tion (12.11) . or its ODE representation (12 .4p . we turn to the main purpose of 
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this paper, which is the characterization of the breakdown of such solutions at 
finite time. By the continuation theorem for autonomous ODEs, if T* < oo 
and [0, T*) is the maximal time of existence of the solution cr, then 

limsup ||cr|| m = oo. (4.1) 

t/T* 

Such a breakdown criterion is not informative enough, as it roughly says that 
"a solution exists as long as it exists". Our main theorem below provides a 
more concise breakdown condition, which is only associated with the stress 
itself, and does not involve any of its derivatives: 

Theorem 4.1 Let cr be a local-in-time solution to ( 12.1 j) in the class 

C([0, T); H m (R 3 )) n C\[0, T); H m -\R 3 )), 

with m > 3. Suppose that [0,T*) is the maximal time of existence, with 
T* < oo, then 

lim / \\<t(-, s)\\l°° ds = oo. (4.2) 

t/T* Jo 



The proof is similar in essence to the proof of the Beale-Kato-Majda 
theorem for the Euler equations j6]. All is needed is an a-priori estimate of 
the form, 



\\<r(;t)\\m<C[t 



f*,jf \\tr(;s)\\ L ~dsY (4.3) 



where C is a continuous function of its arguments, hence (14.11) occurs only if 
(14.21) occurs. The estimate (14.31) is derived in two steps, detailed in the next 
two subsections. 

4.1 A priori estimates for the H m (R 5 ) norm 

Our estimates rely on the following version of the Gagliardo-Nirenberg in- 
equality, 

\\D k f\\L« < CII/II^H^/II^, 

where 1 < k < m and 1 < p < oo, with p = — . With this, we prove the 
following lemma: 
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Lemma 4.2 For f,g,h G H m (M. 3 ) the following triple-product inequality 
holds, 

[ \D^h\\D^f\\D^g\dx<C\\h\\ la 0^^ 

(4.4) 

where (3 < a. 

Proof: We start with the triple product inequality 

/ \D a h\\D^f\\D a ^g\dx< \\D a h\\ \\D^f\\ Lq \\D a ^g\\ LP , 

</R 3 

where 1/p+l/g =1/2. We then use twice the Gagliardo-Nirenberg inequality, 

ii^/iu,<cimi^ 

and 

II n a -Pn\\ r < r , iinii 1_|a ~ /3|/|a| n n a n\\^ a ~^ / ^ 
\\l> g\\LP S \\-u g\\ L \ a _ PMa[ . 

If we choose 

2\a\ , 2\a\ 



q = ——- and p 



\f3\ * \a-P\> 

then (14.4p is obtained. ■ 

We are going to make an extensive use of inequality (14.41) for the case 
where f = er, g = Vu and h = <r (more precisely, /, g, h are components of 
these tensors). Then, combined with the CZ inequality (lA.lljl . we get 

(^(^[^"(Vtt)]) < C'|| CT ||f a| || CT ||^ l/H ||V^||i-J Q - /3|/H , 
which combined with Young's inequality finally gives, 

(D a <r, (D 8 a)[D a -^Vu)}) < C||o-||f a| (||<r|U=o + \\Vu\\ L oo) . (4.5) 

For every \oc\ < m the L 2 -norm of the a-th derivative of cr satisfies the 
energy equation 

^ll^llo + = " (D a (T, D a {u k ■ d k v)) 

+ 2 (D a a, D a [(a(Vu)]) + ^ (D a (T, D a (Vu)) . 

(4.6) 
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The last term is easily estimated using the Cauchy-Schwarz inequality and 
the CZ inequality flATTll . 

(D a a,D a (Vu)) < C\\D a a\\l. 

The middle term can be written as 

(D a cr, D a [(*(Vu)]) = ( Da(T i {D P (r){D a -PVu)) , 

which is a finite sum of terms, each of which can be bounded using (14. 5p . 

Remains the advection term. Because u is incompressible, the term 
Ukdk{D a cr) vanishes, which means that u is differentiated at least once, and 
we can use (14. 5 \ once again. Thus, we obtain the inequality, 

1C W«.|| a 4. 1 

2dt llD CT|lo+ y 

and summing up over all |a| < m, 



~\\D a <r\\ 2 + yII^o-Ho <C[1+ \\<t\\ l ~ + || Vu|U-] \\D a tr\\l, 



dun 1 n n ^ ^ t-i n n ,, i 

TF m + T CT U <C(1+ \\Cr\\ L oo + Vtt 

at A 



A simple integration yields, 



< 



exp 



CI (1 + ||<t||l°° + || Vulvae) ds 



o 



|co||m- (4.7) 



A comment: the energy inequality (14.61) is only formal since we have 
not shown that the H m (M. 3 ) norm of er was differentiate. To rectify this 
delicacy, one has to carry all estimates with the mollified solutions er e , which 
are differentiate in all Sobolev spaces, and take the limit e — > 0, only once 
we have obtained a final estimate for ||cr e || m in the integrated form . 

We will need one more estimate. For all indices i,j, k we have 







dk<Jij+-d k a ij = -d k (uidi<Jij)+dk (au(diUj) + (diUi)aij)+—d k (diUj + d. 



■jUi 



dt A A 

Multiplying by (d k <Jij) 3 (with summation over all indexes) and integrating 
over M 3 we get 



ld..„ ll4 



-/ (dk<Jij) 3 d k (uidi<7ij) dx 
2/ (dkaijf d k [au(diUj)} dx + / {d k ayf d k {d i u j ) dx. 

JR3 A J R 3 
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Using for the first two terms the triple product inequality, and the Holder 
inequality for the third, we get 



1 d 



1 



- — HVer||2 4 + tIIV<7||*4 < C\ || Vcr 



AdV 



X 



L 4 



II 4 



|Vu| 



I cr || loo || VVu ||^4 || V<x| 



+ ||Vcr||3 4 ||VVn|| L 4], 
and after applying once again the CZ inequality flA.llj) 



L 4 



— 1| Vcr|| L 4 < C [1 + ||Vw||loc + ||cr|| L oo] || Vcr|| L 4, 



from which we get 



I V<t|| L 4 < exp 



c [ (l + H 

Jo 



IVtil 



ds 



IVct 



o||l 4 - 



(4.8) 



4.2 L°° estimate of Viz 



So far we have shown in (14. 7[) that if the solution cr breaks down at time T* 
then 



lim / (1 + ||<t|| r oo + II Villi r°o ) ds = oo. 
To complete the proof of Theorem 14.11 it is sufficient to show that 



cr \\l° 



|Vu||loo) ds <C (t, J 



\ct\\loc, ds I . 



Thus, we need to estimate ||Vu||z,°o in terms of || f]| z,°° - Note that the CZ 
inequality provides a bound for ||Vii||lp in terms of ||c||lp for all p (such 
a bound exists within any of the Sobolev W m,p norms), but the prefactor 
is linear in p, hence p cannot be taken to be infinite. Instead, one has to 
perform a more delicate analysis. 

Consider the integral relation f)2.6p between cr and Vu. We split the 
domain of integration into an "outer domain" \y\ > R, a "middle annulus" 
e < \y\ < R, and an "inner disc", \y\ < e, namely, 



Vu 



1 

5za 



tr cr 



1 

lixv s 
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where 



h(x)= [ M^(y):a(x-y)dy 

jR<\y\ 



h(x) = 



M {2 \y) : a{x-y) dy 



e<\y\<R 



J 3 (x) = (P.V.) / M^(y):a(x-y)dy. 

J\y\<e 

Recall that is homogeneous of degree —3 and averages to zero on the 

unit sphere. 

The "outer" integral is estimated using the Cauchy-Schwarz inequality, 



\h(x)\ < 



C 



°" 0- 



The "middle" integral is estimated by taking out the infinity norm of the 

stress 

|/ 2 (a:)| <C|H| L =olog-, 

e 

For the "inner" integral, we exploit the fact that averages to zero on 

the unit sphere to subtract 



so that, 



(P.V.) / MV>(y):<T(x)dy = 0, 

Ay\<e 

h = (P. V. ) / M^\y) : [a{x - y) - <t{x)\ dy. 

J\y\<e 



By the mean- value theorem (this assumes that <x e C 1 (M 3 ), which is the case 
since m > 5/2), we obtain 



\h(x) < j 


f 




\y\<£ 



M p) (y):[(»-M)] 



dy<C 



|y|<£ 



\MW(y)y\» 



i/p 



IVo-l 



Li 



where we used Holder's inequality, and £ represents an intermediate point. 
Since is homogeneous of degree —3, this L p -norm is finite provided 

that 2 — 2p > —1, i.e., p < 3/2, and consequently q > 3. Setting q = 4 and 
combining all three contributions, we get 



(P.V.) / M^\y):a{x-y)dy 

JR 3 



< c 



R 3/2 



|<t|| + || V<t|| L 4£ 



1/4 



|<t||l<x, log 
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It remains to choose R, e such to minimize the bound. Taking 



R 



2IFII0 



and e 



4||<x| 



1/4 



l^lli 00 J Vll^°"lli 4 
we finally obtain 

||Vit||i«, < C ||cr|| LO o (l + log + || || o + log+ || Vcr|| L 4) 

where log + x = max(logx, 0). 

The estimate (14.71) with m = gives, 



(4.9) 



II 47 II o < ex P 
and similarly, from (14.8 



C 



[1 + ||<X|| L 00 + ||Vtl| 



L°° , 



C (1+ IklUoo + ||Vu| 



ds 



ds 



k^ollo- 



I Vcr 



0\\L 4 - 



1 + / (1 + ||0-||i<» + ||Vtt| 



ds 



|| Vcr|| L 4 < exp 
Substituting into (14. 9 p we get 

||V«||ic» < C||<t||£, 
If we define 

M(t) = [ [1 + ||fr||i,oc + ||Vti||ioo] ds and N(t) = / \\cr\\ L oc ds, 
Jo Jo 

then we have an integral inequality of the form 

M'(t) < CN'(t)[l + M(t)}, 
which we readily integrate, 

J [1 + ||o"||l°° + ||Vm||loo] ds < exp J ||o"||x,do ds^j . 

Combining with (14.71) we have proved Theorem 14.11 

Note that (14.71) together with (14. 10p yield a doubly-exponential bound, 



(4.10) 



< exp 



C exp ( C I ||<r||i°° ds 



In this context, it is noteworthy that a similar doubly-exponential bound 
was derived in BKM for the Euler equation. Ponce in [17] derives a singly- 
exponential bound using as control parameter the L°° norm of the deforma- 
tion tensor, rather than the vorticity as in the original BKM paper. 
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5 Discussion 

In this paper we derived a breakdown condition for solutions of the Oldroyd-B 
model in IR 3 in the limit of zero Reynolds number. This condition is analogous 
to the breakdown condition of Beale-Kato-Majda for Newtonian fluids. It is 
noteworthy that the elastic relaxation time, A (which in non-dimensional 
formulations is the Weissenberg number) plays no role in our analysis. In 
fact, nothing changes if we set A = oo, or, alternatively, set A = oo but retain 
the ratio u p X constant (which corresponds to the Kelvin- Voigt model for a 
viscoelastic solid). 

The main implication of our result is that the efforts toward a global-in- 
time well-posedness theory should focus on the control of the L°° norm of 
the stress. 

In a more physically realistic setting, one should consider the same prob- 
lem in a bounded domain Q. Let us assume for simplicity homogeneous 
Dirichlet boundary conditions for it, and a sufficiently smooth boundary dQ. 
The Calderoh-Zygmund inequality holds in this case (see e.g. [TSJ US])) so 
that local-in-time existence can be proved as for the infinite domain. Differ- 
ences arise in the proof of the BKM criterion, where we need an L°° estimate 
for Vit. For a bounded domain, the integral representation (12. 6p can be 
rewritten with the same kernel integrated over Q, plus a boundary-term con- 
tribution, exactly as in the Poisson representation formula (see [IS]). Another 
alternative is to use the following Green representation formula, 



where Qq is a singular kernel that depends on f2. The existence of Qq along 
with pointwise estimates has to be proved using elliptic regularity theory (see 
Ferrari [19j for details). 

A Some inequalities 

In this appendix we list a number of inequalities used repeatedly in Sections [3] 
and HI We recall that H m (M. 3 ) denotes the Sobolev spaces of scalar, vector 
and multi-dimensional tensor fields, with the corresponding norm || • || m . The 
L 2 -norm, which coincides with the if - norm is denoted by || ■ ||q. The L°° 
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norm is denoted by || • Weak derivatives are denoted by D a , where 

a = («i, a2, 0:3) is a multi-index. 

Mollifiers For / G L P (R 3 ), 1 < p < 00, the mollification operator, J £ , is 
defined by 

(J*f)(x) = hf <t>(^-)f(y)dy, (A.i) 

£ JR3 \ 6 J 

where : 1R 3 — > R is a radially symmetric, positive, compactly supported 
C°° function, satisfying J R3 (j)(x)dx = 1. The important properties satisfied 
by J £ are 

• J e commutes with (distributional) derivatives. 

• J £ is symmetric with respect to the L 2 (M 3 ) inner product. 

• J £ : H m (R 3 ) -> iJ m (R 3 ) n C°°(M 3 ). 

• There exists a C > such that for every / G H m (R 3 ), 

||J £ /-/IL-i<C £ ||/|U. (A.2) 

• For every / G H m (M. 3 ) and k > 0, 

II Je/IUfc < ^ll/IU- (A.3) 

• For every / G if m (R 3 ) and multi-index a, 

ll^/IU- < ^ ll/llo. (A.4) 

Banach algebra property of Sobolev spaces For m > 3/2 the Sobolev 
space i7 m (IR 3 ) is a Banach algebra, i.e., there exists a constant C > such 
that for all f,g£ H m (R 3 ), 

||/^IU<C||/|U|b|U. (A.5) 
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Sobolev calculus inequalities For all m G {0, 1, 2, ..} there exists a con- 
stant C > such that for all f,g G L°°(M 3 ) n H m (R 3 ), 

ll/^IU < C (ll/IU~II^IU + ll/IUII&IU~) - (A.6) 

For m > 3/2 the Sobolev embedding H m (R 3 ) C L°°(IR 3 ) implies that this 
inequality is valid for all /, g G H m (R 3 ). Moreover, for / G H m (R 3 ), g G 
ff m - 1 (R n ), and |a| < m, 

|p a (/#) - /D a p||o < C [|| V/|U« \\g\\ m ^ + \\f\\ m \\g\M • (A.7) 

This inequality holds even if each of the two subtracted terms on the left 
hand side is not in L 2 (M 3 ) [BJ. 

Gagliardo-Nirenberg inequalities The classical Gagliardo-Nirenberg in- 
equality is 

IPVllw < C\\f\\l7 e \\D™f\\ s LP , 
where 9 = k/m G (0, 1), 1 < p, r < oo, and 

1 _ 6 1-6 

q p r 

For the particular case r = oo we have 

< C\\f\\)£ /m \\D m f\\ k £ m , (A.8) 
where 1 < p < oo and 1 < m < k. 

Sobolev embedding theorems ( |20j Ch. 5, [21] p. 168) Let m > 1 

an integer and 1 < p < oo. Then: 

If i _ a > o then W m ' p (R n ) C L 9 (M n ) with i = I - a 
If I _ a = o then W m -P(R n ) c L q (R n ) , Vg G [p, +oo) 
If I _ s < o then ^ m,p (M n ) C L°°(M n ) 

In this paper we use extensively the third embedding for p = 2 and n = 3, 
i.e. 

ff m (IR 3 ) C L°°(R 3 ) if m > -. 
In fact, this is a continuous embedding, i.e., 

||/IU~<C||/|| m . (A.9) 
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The C alder on- Zygmund inequality Let K : M 3 — > M. be a homogeneous 
function of degree —3 that averages to zero on the unit sphere, and for 
/ G L q (R 3 ), q>2, define 

g(x) = (P.V.) [ K(y)f(x-y)dy. 

Then g G L q (M. 3 ) and 

lb||M<Cg||/|| w . (A.10) 

Because the constant grows unbounded with q, this inequality does not carry- 
to the L°°-norm; see, for example, Stein [22] for a proof. In particular, since 
the same kernel relates D a f and D a g, it follows that / G H m (M. 3 ) implies 
g G # m (M 3 ) with 

\\g\\m<c\\f\\ m . (A.ii) 
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